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Abstract 

I consider the partition function of the inhomogeneous 6-vertex 
model defined on the n by n square lattice. This function depends 
on 2n spectral parameters Xi and yi attached to the horizontal and 
vertical lines respectively. In the case of domain wall boundary con- 
ditions it is given by Izergin-Korepin determinant. For q being a root 
of unity the partition function satisfies to a special linear functional 
equation. This equation is particularly good when the crossing pa- 
rameter rj = 2tt/3. In this case it can be used for solving some of 
the problems related to the enumeration of alternating sign matrices. 
In particular, it is possible to reproduce the refined ASM distribution 
discovered by Mills, Robbins and Rumsey and proved by Zeilberger. 
Further, it is well known that the partition function is symmetric in 
the {x} and as well in the {y} variables. I have found that in the 
case of r] = 27r/3, the partition function is symmetric in the union 
{x} U {y}l This nice symmetry is used to find some relations between 
the numbers of such alternating sign matrices of order n whose two 
'1' are located in fixed positions on the boundary of the matrices. Fi- 
nally I derive the equation giving 'top-bottom double refined' ASM 
distribution. 



1 



1. Six- vertex model with domain wall boundary and 
Izergin-Korepin determinant 

Let us consider the inhomogeneous six-vertex model on a square lattice. 
The states of the model are given by assigning arrows to each edge of an 
n X n square lattice so that at each vertex, two arrows go in and two go out 
{ice condition ). Spectral parameters {xi, X2, and {yi,y2, ■■■,yn} are 
attached to the horizontal and vertical lines respectively. Edges point inward 
at the sides and outward at the top and bottom as in figure 1. With these 
boundary conditions which are called domain wall boundary conditions 
the states of the model are in in one-to-one correspondence with a certain 
set of matrices, called alternating sign matrices . 

Xi — , ... . — 

X2 — , . . . , 



yi z/2 yn 

Figure 1 

There are six possible configurations of the arrows on the edges for a given 
vertex as in figure 2. 



= a 



= c 



Figure 2 



The Boltzmann weights which are assigned to every vertex of the lattice 
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have the form: 

a[x-y) = : , h{x~y) = , 1 

sm 7] sm 7] 

c{x-y) = 1. 

The different letters in Figure 1 correspond to the different functions intro- 
duced in equation (0). The x and y are the horizontal and vertical spectral 
parameters, respectively, which depend on the vertex position. The t], some- 
times called crossing parameter has the same value for all vertices. 

The weight of a state of the model is the product of the weights of its 
vertices, and the partition function (state sum) is the total weight of all 
states. 

Let Zn{xi,X2, ■■■,Xn',yi,y2, ■■■,yn) be the resulting state sum. Izergin j2] 
(using the above mentioned work of Korepin ) found a determinant rep- 
resentation for Z. 

Zn{{x},{y})= Yl sin(xi +?7/2) ]J sm{x^ - yj - r]/2) x 

^ detM({x},{|/}) 

ni<j<«'<n sin(xi — Xii) ni<j<j'<n sin(?/j — yji) 

where the entries of an n x n matrix M{{x}, {y}) are 

1 

sin(xi - yj + rj/2) sin(xj - yj - r]/2) 



(2) 



^M = — 7Z , z TT^' z,j = l,2,...,n. (3) 



Let us list some properties of the state sum which are used in this paper 
[see, for example, ^ or 

The function Zn{{x}, {y}): 

• is symmetric in the {x} and in the {y}. 

• is a trigonometric polynomial of degree at most n — 1 in m, with vr being 
quasiperiod: 

where u is one of the variables xi, X2, yi, 7/2, ■■■,yn- 



3 



• is left unchanged by transformation 

x^7r/2-x, y -y, 77 ^ vr - 77, 
while function b (see equations dH)) changes its sign. 

2. The basic equation for the state sum 

Let us fix the crossing parameter 

271 

^= iV' 

where is an integer, > 3. 

Consider the set of the matrices M^^^ given by equation (jSj) with the 
spectral parameter xi shifted by kr]: 

M^'') =M{xi + kr],x2,...Xn,yuy2,...yn), k = 0,...,N-l. (4) 
The variable xi occurs only in the upper row of the matrices M^^^: 

M^^^ = - (5) 

^'^ sin(xi - yj + {k + 1/2)77) sin(xi - yj + (A; - 1/2)7/)' ^ ^ 

j = 1,2, ...,n. 

Consider the sums of entries of the upper row over the set of M^^) 

Af-l 

Y.Mi'^l J = 1,2,. ..,7.. (6) 

A;=0 

Inserting (jH]) into ^ and using a simple identity 

sin 7/ expi{u — r]/2) expi{u + ri/2) 

sm{u + T] / 2) sm{u ~ T] / 2) sm{u — 'rj/2) sm{u + ri/2) ' 

we obtain that all sums (jHI) are equal to 0. Since detM^'^) for every j linearly 

(k) 

depends on Mi J this implies that 



detM^^^ = ^ detM(xi + kr], X2, yi, 1/2, -l/n) = 0. (7) 



AT-l N-l 
k=0 k=0 
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Using representation (j21) for the state sum, we immediately obtain the basic 
equation 



\ Z{xi + kr], X2, ...Xn] yi, 1/2, ■■■yn) H sin(xi -Xi + kri) 

k=0 ^ 1=2 



k=0 
n 

n 



sin(a;i — yj + {k + l/2)ri) sin(xi — yj + {k — l/2)ri) 



-1 



This formula appears to be useful for the study of the partition function Z. 
In this paper we limit ourselves to the case 



N = 3, r] = 27r/3, 



Using an identity 



sm[u + —) sm{u — —) smu 



sin 3m, 



we obtain a much more handy equation 



2 r 

J2\Z{xi + 2nk/3,X2, ...Xn]yi,y2, ■■■yn) 

k=0 ^ 

n n 

Y[ sin(a;i — Xi + 27rfc/3) Y\ sin(a;i — yj + 2nk/3) [■ = 0. 

i=2 3=1 



Let us write 

u for xi, 

Ui for Xj+i, i = 1, n — 1, 
and Ui+n-i for yi, i = 1, ...,n. 

After introducing the function / defined by 

2n-l 

f{u) = Z{u) n sin(M-Mi), 



(8) 



(9) 



(10) 



we can rewrite equation (jH)) in the following form: 

Ott A-TT 

f(u) + fiu + -) + f{u + -) = 0, 



(11) 



where, for brevity, we write f{u) and Z{u) rather than f{u,Ui, .■■U2n~i) and 
Z{u, Ml, ■■■U2n~i) suppressing the variables Ui, i = 1, ...,2n — 1 which remain 
unchanged throughout. 

The point is that equations (jlOj) and can be used to find 

the state sum Z{u) up to an arbitrary constant. 

Using the properties of the state sum hsted in Section 1 we obtain that 
the function f{u) is a trigonometric polynomial of degree at most 3n — 2 
which satisfies: 

f{u + n) = {~irf{u). (12) 
This polynomial is expressed by a finite Fourier sum: 

3n-l 

fiu) = E ^fc exp{z(3n - 2k)u}. (13) 

k=l 

where coefficients bk, k = l,2,...,3n — 1 depend on the variables Ui, i = 
l,...,2n — 1. Inserting this expression into (|11|) we find that every third 
coefficient is zero: 

h^ = 0, K = l,2,....n-1. (14) 

Let Ui 7^ Uj for i j. In this case the remaining 2n coefficients can 
be found up to normalization by solving a homogeneous system of linear 
equations^: 

3n-l 

fiuj)= exp{i{3n - 2k)uj} h = 0, j = 1, 2, ...2n - 1, (15) 

k=l,k^3K 

since the rank of the coefficient matrix is 2n — 1 for generic values of Uj . 

The resulting f{u) is given (up to an arbitrary constant) by the determi- 
nant 



P = det 


^3n-2 
t" 


.3n-2 
i^l 

j3n-4 
f^l 

V-3n-8 
f^l 


.3n-2 
'^2 

+3n-i 
l2 

+3n-8 
l2 


.3n-2 
■ ^2n-l 
7-371-4 

• ''2n-l 
7-371-8 

• ''271-1 




^2-3n 


.2-3n 
1^1 


,2-3n 
'^2 


.2-3n 
• ^2n-l 



This problem is a special case of interpolation problem. 



6 



where t = exp{i u) and tj = exp{iuj), j = l,2,...,2n — 1. Indeed, let us 
evaluate the determinant by the left column expansion: 

3n-l 

Pi'^) = exp{i{3n — 2k)u} bk , 

where bk are the corresponding cofactors depending on the variables uj, j = 
1, 2, 2n — 1. We see first that the determinant P{u) is expressed by a finite 
Fourier sum of the same kind as (fT^ . Second, condition (fT^ is satisfied. At 
last it is obvious that P{u = Uj) = 0, j = l,2...,2n — 1 and hence the 
coefficients bk satisfy the system of equations (P3j) . Therefore the ratio 

P{u)/f{u) = C{Ui, U2, U2n-l) (16) 

does not depend on the variable u. 

3. The symmetry of the state sum in the spectral parameters 

Inserting definition ()10|) into equation ()16|) we obtain 

P{u;Ui,U2, ...,M2n-l) = C{Ui,U2, ...,U2n-l) 
2n-l 

Y[sm{u-Uj) Y[ sm{uj - Uj) Z{u;ui,U2, ..■,U2n-i)- (17) 

j=l l<j<i'<2n-l 

Picking out the factor ni<j<j'<2n-i sm{uj — u'j) we emphasize the symmetry 
of the functions entering equation p7|) . Let us return to the parameters {x} 
and {y} (see definition Q). Both the left hand side of equation (fTTj) and 
the product on the right hand side are antisymmetric in the union {x} U {y}. 
Combining this property with the above-mentioned symmetry of the state 
sum Z({x},{y}) in the {x} and in the {y} separately, we obtain that the 
unknown function C has to be symmetric in the {x} and in the {y} as well. 
From the other hand this function does not depend on the variable u = xi. 
Consequently it does not depend on any of the the variables Xj, i = 1, 2, n. 

Returning to the beginning of Section 2 and using the left column instead 
of the upper row, we obtain that the function C does not depend on variables 
yi, i = 1,2, ...,n as well. Hence C is a constant which possibly depends on 
n. So we get the nice result that 

for 7] = 27r/3, the partition function Z{{x},{y}) is symmetric in 
the union {x} U {y}. 
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4. The refined ASM conjecture 



Each of the further sections address some aspects of the ahernating sign 
matrices (ASMs) of size n, which are in one-to-one correspondence with the 
states of the model we consider ini^El- The number of ASMs, for example, 
is equal to the state sum Z with a = 6 = c = 1. 

Consider the case when all spectral parameters, except one are equal to 
zero: 

X\ = u, 

Xi = 0, i = 2,3,...,n, (18) 
Z/i = 0, j = l,2, ...,n. 

In the upper row of the lattice the Boltzmann weights are given by: 

2 2 
a(M) = sin(7r/3 + m), b{u) = sin^ir / 3 — u) , (19) 
V 3 V 3 

c{u) = 1. 

All the Boltzmann weights in the remaining rows are equal to 1. 

Let us divide the set of the states into n classes according to the column 
number 'r' of a (unique) vertex of the upper row where the horizontal arrows 
change their direction. Let A{n, r) be the number of the states of the r^^ 
class, then the partition function (the state sum) is given by: 

n 

Z{u) = Mn, r) a'-'iu) 6"-"(m). (20) 

r=l 

In terms of ASMs, it is also the number of ASMs of order n whose sole '1' 
of the upper row, is at the r*^ column. The wonderful story of these numbers 
can be found in the book by Bressoud [7j. Mills, Robbins and Rumsey 
discovered a nice formula for A{n,r) jHlEl- Their conjecture which is known 
as the refined ASM conjecture was proven by Zeilberger [Hj, who found the 
state sum Z{u) extending Kuperberg's method pij. They both began with 
Izergin-Korepin determinant and solved a difficult problem, since both the 
determinant and the product in the denominator of equation 'badly' 
vanish. 

We also intend to find A{n,r) using equations (fTUI) . (fTT|) and general 
properties of the state sum Z{u) listed in Section 1. It does not take a 
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lot of time, because the corresponding problem has been already solved in 
paper [TU] . 

Let us return to equation ((TU). We have estabhshed that the function 
f{u) can be presented by the finite Fourier sum ()13|) . Taking into account 
that some of the coefficients bk are equal zero (see equation we can 

rewrite it to the form 

n-l 

f(u) = J2 (6+e'(4-3n+6m)n ^ ^- g-^(4-3n+6m)^.^ ^ 
m=0 

In contract with the previous case we cannot use equations (fT^. Instead we 
use the fact that in the last case (see equation (HH)) the monomial sin^"~^ u 
divides this function due to definition ((TU)): 

/(m) = Z(u)sin2"-iM. (21) 

This condition fixes the coefficients 6^, m = 0, 1, n — 1 up to an arbitrary 
constant. Indeed, the first 2n — 2 derivatives of the function f{u) have to be 
zero at M = 0: 

/«(«)k=o = E + (-1)'&™)(4 - 3n + 6my = 0, 

m=0 

/ = 0, l,...,2n-2. 
If / is even then we obtain the system 

E + &„.)(4 - 3n + 6mf^ = 0, A = 0, 1, n - 1, 

m=0 

which (from the non-singularity of the Vandermond matrix) implies 

^m + ^m = 0, m = 0, l,...,n- 1. 
The odd derivatives give 

n-1 

E &m(4 - 3n + 6m)2^+^ =0, A = 0, 1, n - 2. 
This system is equivalent to the condition that the relation 

n-1 

bi{4: - 3n + 6m) p{{4: - 3n + Gmf) = (22) 

m=0 
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is valid for all polynomials p{x) of degree n — 2. Consider the n — 1 special 
polynomial of degree n — 2: 

n-l 

Pf^i^) = n (a; - (4 - 3n + 6m)^), /i = 1, 2, n - 1. 

m=l,m^fi 

Inserting these polynomials into (j^^ we find simple relations connecting 6+ 
with Bq. It is possible to write the answer in terms of binomial coefficients, 
namely, 

n-l / 4 \ / 2 

n — ^ \ I n — 



m=0 



/H= const Y.[ ^MU_^ii j sin(4 - 3n + 6mK (23) 



This function satisfies the second order linear differential equation 

/" - Q{n - 1) cot3n /' - (3n - 2)(3n - 4)/ = 0. (24) 
We return now to equation ((201) • Combining it with equation (^1]) we obtain 

f{u) = sin^"-^ u sin'^-i(M + 7r/3) ^(^, ^) (^)" (25) 

r=l \a('Uj/ 

Inserting the last equation into ODE (j24|) and changing variable to 

^ h{u) sin(7r/3 — u) 
a{u) sin(7r/3 + m) ' 

we obtain the second order linear differential equation 

t(l - t)A" + 2(1 - n - t)^ + n{n - l)A = 0, (26) 
for the generating function: 

n 

A{t) = J2Ain,r)e-^. (27) 

r=l 

Inserting this sum into differential equation ()26|) we get a recursion relation: 

(2n — r — 1) r A(?t,, r + 1) = {n — r) {n + r — 1) A{n, r), 
which is satisfied by 

_ A(n, 1) (n + r-2)!(2n-r-l)! 

" (2n-2)! (r-l)!(n-r)! ' ^^^^ 
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It is clear that the number of ASMs of order n with '1' at the left upper 
corner A{n, 1) is equal to the total number of ASMs of order n — 1: 

A{n,l) = An-i. 

Using this equality and equation (j28p we obtain a recursion relation for the 
total number of ASMs: 



r=l 



{2n-2)\ fr{ (r- 1)! {n - r)\ 



The sum on the right hand side of this equation can be found by comparing 
the power series expansions on both sides of the next identity ^ 

a;^-" (1 - a;)~2" = (1 - a;)"" x x^'" (1 - x^. 

As a result we get the famous recursion: 

An+i (3n + l)!n! 



An (2n+ 1)! (2n)!' 

Using y4i = 1 we find the total number ofnxn ASMs, which was conjectured 
by Mills, Robbins, and Rumsey El and then was proved by Zeilberger JT] 
and Kuperberg jlj. 

5. Double ASM distributions 

Consider now a more general case of two nonzero spectral parameters: 

xi = u and x„ = it, (29) 

which are are attached to the upper row and the lower row, respectively. 
Here the analogue of equation for the state sum becomes 

n 

Z{n- u,u)= B[n- r, f ) a'-^[u) h'^-^u) ¥~\u) a^'-^u), (30) 

r,f =1 

where B{n, r, f) is the number of ASMs of order n whose sole '1' of the upper 
row, is at the r*'^ column and whose sole '1' of the lower row, is at the r^^ 
column. 



^It is also possible to have a look at standart textbooks on combinatorics. 
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According to the symmetry, announced in Section 3, we have to get the 
same value of the state sum by fixing the non-zero spectral parameters as 
follows: 

xi = u and yi = u. 

In this case the non-zero parameters are attached to the upper row and the 
left column and we can present the state sum in a slightly different way 

Z(n; u, u) = A{n; 1) + a{u - u) 

Cin; r, f) a''-\u) a^-\-u) (31) 

r,f=2 

where C{n; r,r), r, f > 2 is the number of ASMs of order n whose sole '1' in 
the upper row, is at the r*'^ column and whose sole '1' in the left column, is 
at the f*'* row. A{n, 1) is the number of ASM with '1' at the left and upper 
corner. 

Equating (jHUj) and (jSH), using equalities a(±£t) = b{^u) (see flTIJ|) ). and 
dividing the resulting equation by (&(m) a{u))^~^, we obtain 

+ 6(^)a(«),3^,^^^'"'"H6(«)J \a{u)) ' ^^^^ 

Now we have to use the identity: 

a{u — u) a{u) ^ h{u) ^ 
h{u) a{u) h{u) a{u) 

Writing t and t for the ratios a{u)/h{u) and h{u)/a{u) respectively we can 
present equation (j32p in more transparent way 

n 

J2 B{n- r, f ) f-^ f-' = A{n, 1) + 

r,r=l 

n 

+ {t + i-l) ^ C(n; r, f ) f 

r,f=2 
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Equating the coefficients we get 



A{n,l) =C{n;2,2), r = r = 1, 

B{n; r, f) = C(n; r, f + 1) + C(n; r + 1, f) — C(n; r + 1, f + 1), 
r + f > 2. 

The first equation is a consequence of a simple bijection [T2I. The second 
equation for f = 1 and r > 1 gives the recursion relation for C{n; r, 2): 

C(n; r, 2) - r + 1, 2) = r, 1) = v4(n - 1, r - 1), 

where y4(n — l,r — 1) is the number of ASM matrices (n — 1) x (n — 1) 
for which the (unique) '1' of the ffist row is at the (r — 1)*^ column. This 
recursion relation can be obtained by a simple bijection as well. I reckon that 
the remaining relations (r, f > 1) are new. 

In the final part of this section we intend to reduce the double ASM 
distribution B{n;r,r) to the numbers A{n,r) (see equation ()28j) . In the 
case of two nonzero spectral parameters (see equation ()29|l ) we can rewrite 
equation (fTUj) as: 

fdoubieiu) = Z{n; u, u) sin^""^ u sm{u - u). 

This function is the trigonometric polynomial of degree at most Sri — 2. It 
can be found up to an arbitrary constant by using equations (fTTj) and (fT^ . 
Let us consider the combination 

fdoubieiu) = const{f{u)f'{u) - f'{u)f{u)), 

where the function f{u) was found in Section 4. This combination satisfies all 
above-mentioned conditions and hence solves the problem. Inserting ()25p into 
the last equation we obtain after tedious but straightforward calculations: 

A r~f 1 H(t)G(i) - H(i)G(t) , , 

^ B{n; r, r)t"-'^f-^ = const ^ ' ^ ' (33) 

r,r=l t — t 

We use two new functions: 
G(t) = (l-t)A(t), 

H{t) = (n- 1/2) (1 + t) A{t) + (1 -r2)(l -t)(l/2-t) A{t) - 
-{l-t)it' -t+l)A'{t), 
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where A(t) is the generating function defined by equation (j77|) . Equation 
leads to a nonhnear relation between the 'double refined' ASM distributions 
and the numbers appearing in the refined ASM conjecture (theorem!): 

B{n; r + 1, f + 1) — B{n; r, f) = {A{n — 1, r) [A{n, f + 1) — A{n, r)] + 
A{n -l,r) [A{n, r + 1) - A{n, r)]}/A{n, 1). (34) 
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Note added to the second ArXiv version of this paper 

Combining equation (jl7|) with equation (j2)) we can relate the determinant 
P with the Izergin-Korepin determinant up to a constant multiplier. The 
precise form of this relation can be derived ^2] by using the rj = 27r/3 
specialization of the theorem 4.2 of paper |14j . 
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